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At quantum critical points (QCP) 0, & H 0, H, 0, @] there are quantum fluctuations 
on all length scales, from microscopic to macroscopic lengths, which, remarkably, can 
be observed at finite temperatures, the regime to which all experiments are necessar- 
ily confined. A fundamental question is how high in temperature can the effects of 
quantum criticality persist? That is, can physical observables be described in terms of 
universal scaling functions originating from the QCPs? Here we answer these questions 
by examining exact solutions of models of correlated systems and find that the tem- 
perature can be surprisingly high. As a powerful illustration of quantum criticality, we 
predict that the zero temperature superfluid density, p s (0), and the transition temper- 
ature, T c , of the cuprates are related by T c oc p 3 (0) y , where the exponent y is different 
at the two edges of the superconducting dome, signifying the respective QCPs. This 
relationship can be tested in high quality crystals. 



Do quantum critical points (QCPs) provide a powerful framework for understanding complex correlated many body 
problems? Do they shed new light on quantum mechanics of macroscopic systems, providing, for example, a deeper 
understanding of entanglement? Answers to many such questions require a precise recognition of the experimentally 
observable regime of quantum criticality. An important class of QCPs are analogs of classical critical points, but in 
a dimension higher than the actual spatial dimension of the system; that is, a quantum critical point in d spatial 
dimensions is equivalent to a classical critical point in (d + z), dimensions, where z is the dynamic critical exponent. 
We shall be primarily concerned with 2 = 1, although an example involving z ^ 1 is given below. This extra dimension 
(imaginary time) is the sine qua non of the effect of quantum mechanics. 

It seems paradoxical that often the region of classical critical fluctuations is small and tends to zero as the tem- 
perature, T, tends to zero, while the region of quantum critical fluctuations fans out as T increases. The reason is 
that the quantum critical region is determined by the dominant microscopic energy scale in the Hamiltonian, whereas 
the classical critical region is determined by the ratio of the transition temperature, T c , to the same scale raised to a 
positive power (the dimensionality). Since this ratio is usually small compared to unity (it is 10~ 5 in a conventional 
superonductor), the classical critical region is small as well. As a classical critical point is driven to zero by tuning 
a parameter to reach the QCP, T c itself vanishes as £~ z , provided the T = spatial correlation length, £, is large 
compared to the lattice spacing, further amplifying the contrast between classical and quantum criticalities. 

How high in temperature does the effect of a QCP persist? Consider perhaps the simplest possible example of a 
quantum phase transition described by the Hamiltonian of an Ising model in a transverse field 8] in one dimension, 

H = -5>i(")- Ai5> 3 «Mn + 1) (1) 

n n 

where the usual Pauli matrices, denoted by the its, are situated on a one-dimensional lattice labelled by n, of lattice 
spacing a. Here we have scaled out an overall energy scale denoted by J. When Ai ^> 1, the ground state is that of a 
classical one-dimensional Ising model, which is ordered at zero temperature. The transverse field represented by the 
Pauli matrix a\ introduces quantum fluctuations and tries to restore the broken Z 2 symmetry. In the extreme limit 
Ai — > 0, the system is a collection of isolated two-level systems and the ground state of each is the symmetric linear 
superposition of up and down spin states. It is well known that there is a QCP at the value of the dimensionless 
coupling constant Ai = 1 , belonging to the universality class of the famous two-dimensional classical Ising model 
solved by Onsager. The scale J enters as soon as we attempt to calculate the free energy density at finite temperature 
T, which is (setting the Boltzmann constant ks = 1) 
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where — 1/T, and <I> is the dimensionless scaled free energy. The quantity A^ denotes the excitation spectrum of 
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the Jordan- Wigner fermions 

A k = 2 (1 + \\ + 2\iCoska) 1/2 (3) 

and eo is the ground state energy density. The quantity c = 2Jay\i/h is the velocity of the elementary excitations 
at low energies and is clearly non-universal. 

For the validity of QCP at finite temperatures, the free energy density should satisfy the constraints of finite size 
scaling 9], namely it should be expressible in terms of a universal scaling function $ s . In the neighborhood of the 
QCP, this function has as its argument the ratio of two fundamental lengths and approaches a pure number in the 
limit T — > 0, when tuned to the exact critical point. The extent to which $ approximates $ s defines the quantum 
critical scaling regime; 5> s can be extracted from Eq. 0and satisfies 

/o-eo^-£$ s ([W£oo] lA/ ), (4) 

where L T — hc/T is the thickness in the imaginary time direction and £oo oc 1/|1 — Ail" is the correlation length when 
T = , that is, when the thickness L T = oo; the exponent v = 1 for the present model. 

If we tune to exact quantum criticality, that is, Ai = 1, then it is known that [ToL Hl| <E> S (0) = 7r/12. Since this is 
the best chance for quantum criticality to survive at finite temperatures, we test our free energy at this point. The 
result is shown in Fig. ??. It is quite remarkable that scaling holds to a temperature as great as J/2. A large overall 
energy scale helps ensure quantum criticality at high temperatures. 

There are always irrelevant operators, however, whose effects must disappear before scaling can be observed. The 
above Hamiltonian was constructed to have no irrelevant operators. We now turn our attention to a modification 
of the above Hamiltonian by adding an irrelevant operator involving a three-spin interaction term. This term is 
generated in the first step of a real space renormalization group procedure [12j; here we allow its coefficient to be 
arbitrary. The new Hamiltonian is 

H = - ^o-i(n) [1 + A 2 cr 3 (n - 1)03 (n + 1)] - Ai ^ a 3 (n)a 3 (n +1) (5) 

n 

The only change now is that the transverse field is modified by the nearest- neighbor spins. We have solved this 
Hamiltonian, once again, by Jordan- Wigner transformation. The fermionic spectrum 

A k = 2(l + A^ + A^ + 2Ai(l-A 2 )cosfca-2A 2 cos2fca) 1/2 (6) 

vanishes at ka — ±7r , when A 2 = 1 — Ai , and at k = 0, when A 2 = 1 + Ai, where Ai, A 2 > 0, thus defining two critical 
lines (see Fig. EJ - 

The introduction of the irrelevant operator with the coupling A 2 does not change the universality class because the 
correlation length exponent is still given by v = 1 all along the critical lines. It is sufficient to examine the free energy 
density for A 2 = 1 — Ai and < Ai < 1. It is 
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where he = 2Jay / 4A 2 + Ai(l - A 2 ) = 2Ja\\ 1 - 2| . In Fig.|3]we plot the function $(Ai, A 2 = 1 - A x , J/T). The results 
show that as we move along the critical line starting at (A 2 = 0, Ai = 1) and ending at (A 2 = l,Ai = 0), we can 
dramatically alter the regime of validity of quantum critical scaling at finite temperatures. 

We analyze another model that is important in the study of antiferromagnetism and high-temperature supercon- 
ductors, the quantum (9(iV)-non linear er-model Taking d = 2, we consider the partition function 

Z= I Pai(n 2 -l) e -i^^^ 2l( « ! (8) 



where fi — (0, 1, 2), and n is a ./V-component unit vector. We obtain an exact expression for the free energy in the 
limit N — > 00, as in Ref. The only difference is that we compute the free energy numerically without assuming 
T/hc <g; A and follow the same procedure as before: tune to the QCP, gi c , and look for a breakdown of finite-size 



3 



scaling. As shown in Fig. scaling is valid up to a temperature T* w hcA/8. Using spin wave theory and an 
appropriate lattice regularization we obtain T* m J/2, where J is the antiferromagnetic exchange constant of the 
Heisenberg model, which is of order 1000 K in cuprates. 

Does T* change in the presence of additional irrelevant couplings? We add the simplest irrelevant perturbation that 
permits an exact solution in the N — > oo limit, namely 172 n • n. Using Eq.|5| given in the Methods section, we have 
estimated T* at several critical points (<7i c ,<72c) and found that it changes very little as we move along the critical 
surface. Hence the irrelevant coupling gi seems to be truly innocuous, in contrast to the parameter A2 in our previous 
model. The difference is not so surprising: in the latter case wc added a perturbation that altered the low-energy 
excitation spectrum — not only by renormalizing the velocity, but, more importantly, by introducing a secondary 
minimum. The term we added here, on the other hand, has no effect whatsoever in the limit to, k — > 0, leaving the 
original low-energy spectrum intact. A different approach to the present model has also suggested T* w J/2 p^|. 
Our analysis of the leading irrelevant operator and our consideration of the free energy have strengthened this result. 

In the cuprate superconductors, quantum criticality leads to a powerful and testable prediction relating the transi- 
tion temperature, T c , and the T = superfluid density, p s (0). Consider first the case of a pure system with (d + z)- 
dimensional QCPs at the edges of the superconducting dome at dopings x = x c . In general the universality classes of 
these transitions are of course different. Scaling [jllUg leads to T c oc \x — x c \ ZUd +- and p s (0) oc \x — x c \ 2t3d + z ~ rid +* l ' d +* 
for x close to x c (with standard definitions of critical exponents). Therefore, T c oc p s (Q) y , where y is determined by 
the universality class of the transition in question. If (d + z) is such that hyperscaling holds (below the upper critical 
dimension d u ), then y = z/(d+ z — 2), independent of Vd+z- Moreover, if d = 2, even z cancels out from the formula, 
resulting in the superuniversal result y = 1. Remarkably, this is also the Uemura plot, further elaborated by Emery 
and Kivelson If (d + z) > d u , the Gaussian fixed point is stable, and y = zVd+ z /{2/3d+z — Vd+z^d+z) = z/2 

because v d +z = Pd+z = lA Vd+z = 0. 

Since the coupling to nodal fermions is irrelevant [ljj, the QCP in the underdoped regime belongs to the XY 
universality class, with z = 1 and d u — 4. Asymptotically close to the QCP, the transition must correspond to 
d = 3, but there is a regime further away in which the layers are effectively decoupled due to weak coupling between 
them, hence the spatial dimensionality can be set to 2. The details of this crossover depend on various microscopic 
parameters that are beyond the scope of a scaling theory and must be established empirically. Thus, there is a regime 
in which Uemura scaling holds. But the true criticality is given by y = 1/2, modulo weak logarithmic corrections. In 
contrast, the QCP in the overdoped regime reflects a superconductor-to-metal transition, suggesting z > 2 and y > 1. 

Strong disorder and macroscopic inhomogeneity can of course destroy the QCPs and therefore scaling. Of more 
interest to us is weak disorder that modulates x c (analogous to modulation of T c in a classical phase transition) , which 
may not change the scaling arguments. The reason is as follows: a simple extension of the famous Harris criterion 
to the T = quantum problem implies that the exponents of the pure system remain unchanged if Vd+z > 2/d. If 
this inequality is violated, there are two possiblities: (1) the system flows to a new fixed point, in which case the 
scaling arguments remain valid but with new exponents; (2) the transition is seriously rounded and the QCP itself is 
destroyed. As in any problem involving disorder, a priori it is difficult to know which of the two possibilities holds. 
The robustness of the scaling argument is quite striking, however. When [d + z) < d u , Vd+z cancels out, and even z 
cancels out in d = 2; when (d + z) > d u , y is simply z/2. 

Experiments in high quality crystals |l8j do in fact yield a value of y = 0.61, substantially less than unity, in the 
extreme underdoped regime. Moreover, based on our analysis, a plot of T c vs. p s (0) should rise more steeply in the 
underdoped regime than in the overdoped regime; a number of experiments appear to support this assertion at a 
qualitative level 

13 E3 (an effect known in the literature as the "boomerang" effect^ although more data are 
necessary, especially because there are other experiments that do not fit into this picture [22|. 

Clearly the role of irrelevant operators is a limiting factor |23j in observing scaling, but, as the term suggests, it is 
difficult to know ahead of time the set of all such operators. In some cases unexpected subtleties can arise at very low 
temperatures, as in a recent experiment involving LiHoF4 [2^ . In this case, the QCP dictated by the electronic spins 
is forestalled by the hyperfine coupling of electronic and nuclear spins. The true QCP is the one at which both the 
electronic and nuclear spins are entangled. This entanglement shows up as a transfer of intensity from the magnetic 
excitation of electronic origin to soft modes of the coupled electronic and nuclear origin in the 10/ieV range. The 
result is an additional quasielastic peak at the measured temperatures due to thermal decoherence. In systems like the 
cuprates, where the dominant microscopic energy scale (the exchange constant, for example) is very large, such effects 
are not of much consequence, and the QCP determined by the larger scale provides the correct picture. Another 
aspect of quantum criticality is the quantum critical fan 0, 0, 1 1 - ll ] that is determined by the crossover temperatures, 
T x , generically given by T x = A\x — x c \ ZVd +* . The amplitude A controls the extent of the fan, but decreasing the 
product zi>d+z clearly narrows it, hence the scaling regime. 

The extended transverse-field Ising model appears to be an interesting example, where the QCP of higher symmetry 
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crosses over to a QCP of lower symmetry at low temperatures (see the Methods section). Whether such a mechanism 
is effective in a real system like the anomalous normal state of the cuprates is presently not known |25| . Can quantum 
criticality explain linear-in-T resistivity of optimally doped high-T c superconductors to temperatures as high as 1000 
K? It has been argued that this is unlikely |26|. Nonetheless, we have shown that simple considerations of quantum 
criticality at the edges of the superconducting dome can lead to very powerful experimental consequences for cuprates 
that can be precisely tested. 



Methods 



The inverse correlation length of the O(N) non-linear cr-model with an irrelevant operator involving higher 

order gradient 



The inverse correlation length m(gi, g 2 ,t) is determined by the modified saddle point equation 
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where e\ — k 2 + g 2 k A + m 2 , t = T/Hc, and gi = (N — l)gi. The integrands \?i and Vf^ are given by 
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(9) 



Duality transformation 

By a duality transformation |2?| where the dual quantum spins (Ti(n)) satisfying the Pauli spin algebra are located 
at the centers of the bonds of the original lattice, we can rewrite our extended model (Eq. 5) as 

H = -J2 [Am(n) + T 3 (n)T 3 (n + 1) - A 2 r 2 (n)T 2 (n + 1)] , (10) 

n 

which has been recently studied in the context of entanglement j2^. The point (Ai = 0, A 2 = 1) is the unstable isotropic 
XY fixed point; note that (Ai = 0, A 2 = —1) has the same excitation spectrum. For temperatures T > T* oc A (see, 
Figs. |21 and the behavior is controlled by this unstable fixed point. An integrable perturbation induces a flow to 
the Ising fixed point at lower temperatures. In the language of conformal field theory, this crossover is described by 
Zamolodchikov's "c- function" , which takes one from a theory of conformal charge c = 1 to that of charge c = 1/2 29] . 
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FIG. 1: The scaled free energy "^(Ai = 1, J/T) plotted as a function of 2J/T. The asymptote as T — > should be the 
universal number tt/12, if quantum critical scaling holds at finite temperatures. It clearly does so with sufficient accuracy for 
temperatures as large as J/2. 
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FIG. 2: The critical lines of the Ising model in a transverse field with three spin interaction. There is a region of parameter 
space near Ai = and A2 = 1 in which the system will be close to both critical lines, meaning it has low energy excitations 
at both k = and ka = ±vr, with gaps Ao = 2J|Ai — (A2 — 1)| and A w = 2J|Ai + (A2 — 1)| respectively. Ao collapses along 
A2 = 1 + Ai and A T collapses along A2 = 1 — Ai. In either case, the critical exponent is given by v — 1. Under the duality 
transformation described in the Methods section, the multicritical point (Ai = 0, A2 = 1) maps onto the isotropic fixed point 
of the quantum XY model in a transverse field. 
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FIG. 3: The quantity $(Ai,A2 = 1 — Xi,J/T) is plotted as a function of 2J/T for different values of Ai while maintaining 
criticality, that is, A T = 0. From top to bottom, the values of Ai are 0, 0.02, 0.04, 0.06, 0.08, 0.1, 0.2, 0.4, 0.6, 0.8, 1.0. Five 
of these curves are shown in color, with the corresponding excitation spectra displayed in the inset. In every case, the gapless 
excitations near ka — ±7r are responsible for scaling; but as we approach the multicritical point (A2 = 1, Ai = 0), a new set 
of low energy modes appears at k = 0, with the gap Ao/J = 2|Ai — (A2 — 1)| = 4Ai. $ does not approach <& s = 7r/12 until 
T <C Ao, which requires arbitrarily low temperatures as Ao — > 0. Precisely at the point (A2 = 1, Ai = 0), Ao collapses and the 
number of critical modes is doubled, giving rise to a new asymptote <3> s = n/6. 
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FIG. 4: The free energy density minus the ground state energy density, (/ — eo)(ftc) 2 /T 3 , of the O(N) quantum non-linear 
a- model in the limit N — > oo (dotted line). In the scaling limit (T — > 0), this quantity approaches the asymptote —0.153 • • • 
(solid line). Quantum critical scaling holds up to T* « ftcA/8 with 1% accuracy. 



